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The Advantages of Dynamic Stiffness Parameters 
Over Classical Bearing Coefficients 


An ideal model for a rotor/bearing system will enjoy 
several properties. First, it will have parameters that are 
independent of coordinate system orientation. It will be 
capable of handling isotropic as well as anisotropic systems. It 
will encompass fluid as well as rotor dynamic characteristics in 
a single, elegant model. Finally, it will be useful for analyzing 
and predicting all types of stability problems, particularly fluid- 
induced instabilities which are a very common problem in 


rotating machinery. 


This article takes such a model, developed by Don Bently and 
Dr. Agnes Muszynska, and shows how its Dynamic Stiffness 
parameters provide significant benefits over the bearing 
coefficients used in classical models. The relationship between 
Dynamic Stiffness parameters and classical bearing coefficients 
is shown, as well as techniques for transforming between the 


two for a general, anisotropic bearing. 
The Bently/Muszynska Model 


The fluid force model developed by Bently and Muszynska [1] 


Zone dranback of pertains to rotors at low and medium eccentricities. Among 


ae other contributions, this model has clarified understanding of 
such coefficients, 


: : If-exci ibrati h il whirl hip. 
however, is that it self-excited rotor vibrations such as oil whirl and whip. An 


isotropic bearing is approximated by the Bently/Muszynska 


makes prediction of 
fluid-induced 
instability difficult. " 


model as follows: 
-F, = m,(# — j2AQi - Qr) + Dr + (K - jAQD)r 


=0 if the fluid inertia effect m p =0 oris negligible (1 ) 


Here, Fr denotes the fluid force, mr denotes the fluid inertia 
effect, K and D are the fluid stiffness and damping, Q is shaft 


speed, and À is the fluid circumferential average velocity ratio. 
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ELLIPTIC ORBIT AS A COMBINATION OF FORWARD AND REVERSE COMPONENTS 


The rotor lateral response is described by the vector of 
displacement r=x+jy (x and y are horizontal and 


vertical displacements, respectively, and j =~v-1). 


Equation (1) represents a fluid force rotating forward at 
average speed AQ with stiffness K, damping D, and the 
fluid inertia effect mp. Unlike classical bearing 
coefficients, Dynamic Stiffness parameters such as the 
previously mentioned K, D, A (and reverse stiffness and 
damping, K' and D’, to be introduced later) are 


independent of x,y coordinate system orientation. 


"Another drawback of bearing coefficients 
is that they are not independent of the 


coordinate system orientation." 


In general, fluid bearings attached to a rotor are 
anisotropic. Thus, a model which includes the effects of 
asymmetry to analyze a general anisotropic bearing is 
desirable for the most accurate results. A more 
generalized version of the Bently/Muszynska model of 
equation (1) does include anisotropic effects, and is 
shown later in this article. This model, through its use of 
the parameter À, clearly indicates the instability margin, 


even for rotors with anisotropic bearings [2]. 
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FIGURE 1 
Bearing Coefficients 


In practice, some engineers use bearing coefficients [3] 
from manuals or handbooks for rotor dynamic analysis. 
This is particularly commonplace when measured data 
is not available. One drawback of such coefficients, 
however, is that it makes prediction of fluid-induced 
instability difficult. Some have argued that fluid- 
induced instability is not nearly the problem that it was 
fifty years ago. This is not true. While it may show up 
less frequently in operating machines today, that is 
primarily because designers consider it more fully 
in their designs and take measures to avoid it, 
often by compromising other desirable attributes 
of the rotor dynamic system and by using less- 
than-ideal bearing designs in an attempt to 
preclude these instabilities. Because designers 
must be concerned with such instabilities, the 
ability for a model to fully address and predict fluid- 


induced instabilities remains vital. 


Another drawback of bearing coefficients is that they 
are not independent of the coordinate system 
orientation. This can complicate data collection and 
reduction, as excitation and measurement locations do 
not always coincide with the most desirable coordinate 
system orientations. The use of a single complex 


variable, r, instead of bearing coefficients with x and y 


variables, is a better and simpler way to model the 


bearing force and the rotor response. 


Transformation between Dynamic 
Stiffness and Classical Bearing 
Coefficients 


As noted earlier, this article establishes a method of 
transformation between Dynamic Stiffness parameters 
and classical bearing coefficients for a general 
anisotropic bearing. This is developed as a series of 
steps. First, the elliptic response of the rotor due to 
anisotropic bearings is described using the concepts of 
full spectrum, whereby any generalized elliptical 
response can be expressed as a linear summation of 
forward- and reverse-rotating vectors of various 
rotational frequencies. Next, the fluid force in an 
anisotropic bearing is introduced using Dynamic 
Stiffness parameters, followed by the response 
solutions expressed in terms of their forward and 
reverse (i.e., full spectrum) components. Finally, 
Dynamic Stiffness parameters are given based on 
bearing coefficients and vice versa. We then show how 
tluid-induced instability can be readily predicted with 


these Dynamic Stiffness parameters. 
Elliptic Orbit Construction 


With x and y probes installed laterally, as shown in 


Figure 1, any arbitrary elliptical orbit with filtered 
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frequency @ can be constructed in terms of forward 
and reverse rotating vectors. The rotor response with 
frequency @ at the probe locations can be given by 
r= —e’ ejor B oJ ge (2) 
2 2 


where 


Ae” = (ae + jA e) 


1 m 

Be® = a (4e + jA,e’ ‘ 

A,,a, = amplitude (peak-to-peak) and phase lag in Probe x 
with filtered frequency œ 

A,,@, = amplitude (peak-to-peak) and phase lag in Probe y 
with filtered frequency @ 


A,a@ = amplitude and phase lag as shown in Figure 1 
with forward frequency @ 


B,B = amplitude and phase lag as shown in Figure 1 
with reverse frequency -@ 


For synchronous (@=@) and nonsynchronous (@#Q) 
response, the filtered orbit will be elliptic whenever a 
reverse component B exists. As shown in Figure 1, the 
major axis is half of the forward component A and the 
reverse component B ; the minor axis is half of A 
minus B ; and the orbital inclination angle is halt of the 


phase angles a and B. 


Fluid Force in an Anisotropic Bearing 


While an isotropic bearing leads to a circular orbit 


response for a symmetric rotor, an anisotropic bearing 


MAJOR / MINOR STIFFNESS AND DAMPING AXES DUE TO ANISOTROPIC EFFECT 
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FIGURE 2 
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results in an elliptical orbit response. As shown 
previously, the latter causes a combination of both 
forward and reverse responses. To account for this 
effect, additional terms are introduced into equation 
(1). When fluid inertia is neglected, the fluid force in an 
anisotropic bearing, from equation (1), can be written 
as 

-F, = Di +(K-jAQD)r +D'e?°F+K' e? 7 (3) 

anisotropic effect 

where F = x- jy denotes the complex conjugate of r (if 
r rotates at forward frequency @, for example, 
r = peľ), then F will rotate at reverse frequency 
-0, i€., r= ne) ). Reverse stiffness and damping, 
K' and D', cause reverse component B as indicated 
previously. They are oriented at 6, and 6,, and vary 


with the x, y coordinate system orientation (for 


perturbation. Relations between the rotor response as a 
combination of forward and reverse components and 
the perturbation force can be established in terms of 
Dynamic Stiffness parameters. Therefore, these 
Dynamic Stiffness parameters can be obtained through 


direct and quadrature plots versus perturbation speed. 


Transformations between Dynamic 
Stiffness Parameters and Classical 
Bearing Coefficients 


Dynamic Stiffness parameters (K, D, A, me, K, Di, Og, 6p) 
provide a clear description of bearing characteristics, 
especially instability margins. Bearing coefficients, 
however, do not relate the four stiffness components 
Karkar <G Ky) to the four damping components 


XX! XY! YX! 


(Da, Da, D 


XX! xy! Jx! 


Dy, ). Thus, it is necessary to establish 
a transformation between 


these two expressions. 


“It is suggested that the Dynamic Stiffness parameters 


(which give clear insight into instability scenarios), or the 


From equation (3), if the 


bearing coefficients as a function of the Dynamic Stiffness 


parameters, be used in rotor dynamic analysis.” 


example, 6, and 6, would decrease by 20° if the x, y 
coordinate system rotates 10° forward). Major/minor 
stiffness and damping for this anisotropic bearing are 
related to the anisotropic effect, as shown in Figure 2. 
The major stiffness axis is oriented 6,/2 relative to the 
x-axis with a maximum value K + K', and the minor 
stiffness axis 6, /2+90° with a minimum value K- K’ 
The major/minor axes for damping have the same 
properties as those for stiffness. In general, major/minor 
axes for the ellipse of a rotor’s orbital response, 
stiffness, and damping, are neither co-linear nor 
perpendicular. At low speed, however, where damping 
effects are negligible, the major axis for the ellipse of 
rotor response Is the same as, or co-linear to, the minor 
axis for stiffness. Note that these major/minor axes exist 


only if reverse stiffness K’ and reverse damping D' exist. 


Dynamic Stiffness parameters for an anisotropic bearing 


can be determined by performing nonsynchronous 
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x, y coordinate system Is 
used instead of the com- 
plex variable r=x+jy, then bearing stiffness and 
damping matrices can be expressed in terms of 


Dynamic Stiffness parameters as follows: 


E d 

K,, K, (4) 

K-#Q?m, AQD ue cosé, sind, 
-29D  K-HQ*m, sind, —cosd, 


anisotropic effect 


Dy Dy 
D D| 


yx yy 


. 5 
| D eo eof med (5) 


-2AQm, D sind, —cosd, 


anisotropic effect 


Therefore, one can obtain bearing stiffness and 
damping coefficients if the Dynamic Stiffness 
parameters are known [Table 1]. Likewise, one can 


calculate Dynamic Stiffness parameters if the tabulated 


BEARING COEFFICIENTS IN TERMS OF 
DYNAMIC STIFFNESS PARAMETERS 


BEARING COEFFICIENT EXPRESSION 
K (N/m, Ib/in) = K — ¥Q’m, + K' cos ôx 
K,, (N/m, lb/in) = AQD + K' sin 6, 

K,, (N/m, Ib/in) = —AQD + K' sin 6, 
K,, (N/m, Ib/in) = K — VQ'm, — K' cos ôx 
D,.(Nes/m, Ibes/in) = D + D' cos 6, 
D,, (Ne s/m, lbes/in) = 2AQm, + D' sin 6p 
D,, (Ne s/m, lbes/in) = —2AQm, + D' sin 6p 
D,, (Nes/m, lbes/in) = D — D' cos 6, 

TABLE 1 
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bearing coefficients are known [Table 2]. Note 
that bearing coefficients in handbooks are often 
given in nondimensional form. For instance, the 
following nondimensional stiffness and bearing 
coefficients kij and di; are used in a bearing data 
handbook [3]: 


C, C,Q 
k= Ky d; = D; (6) 





where z, 7 stand for index x or y and W, C, , and 
Q denote bearing load (the radial force), 
bearing/journal radial clearance, and rotative 
speed, respectively. Therefore, the fluid average 
circumferential velocity ratio, A, can be 
expressed directly in terms of nondimensional 
bearing coefficients [Table 2]. As is shown in the 
table, when cross damping coefficients Dy and 
D,,, are the same, fluid inertia mpwill be zero. 

This is often the case when 


using tabulated bearing 


DYNAMIC STIFFNESS PARAMETERS IN TERMS OF BEARING COEFFICIENTS 


DYNAMIC STIFFNESS 








coefficients, and indicates a 


negligible amount of fluid 


EXPRESSION : . : 
PARAMETER inertia. Consequently, fluid 
stiffness Ķ will be the average 
i Bi or Bey T Boe when : value of Kx and K . In seals, 
4 = Q(D,,+D,,)> dy + 
z a however, D and D, are no 
l ; ; y Ix 
nondimensional coefficients k., d., are used 
j> Gi : longer the same (often Dy = 
ee 1 (K, ; K,) fi 1 (K, 7 K) Dy =D; : Dig so that fluid inertia “i 
2 4 Dt), : cannot be neglected. Fluid 
lin ay : damping D is always the 
° ° j = pu xx + 
Da sean) 2 ») average value of D, and D,, . 
| i (D, =D (Dj +D As long as Ky, # Ky, or Ky, # 
m, (kg, Ibe s* /in) = —— . -K., reverse stiffness K’ will 
4 K,, — K, yx! 
exist. Likewise, if D,,, # D, or 
K' (N/m,Ib/in) = 2 (Ke — K,,) + (K; + K Do # Dy reverse damping 
D' will be nonzero. Their 
D' (Ne s/m, Ibes/in) = l l ientati 
(N-s/m, lbes/in) > \(D, D) + (D, + D,,) orientation can be calculated 
and is related to the 
ô, (deg) = tan”! Ky + Ay major/minor stiffness or 
K,. — Ky, 
damping axis. 
D +D 
On (deg) = tan! O N o 
Dax — Dy 


TABLE 2 
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Summary 


This article has established several compelling reasons 
for the use of Dynamic Stiffness parameters (K, D, A, me, 
K', D’; 6x, 65) instead of less meaningful classical 
bearing coefficients K og Ky, Ky 1 Ky 1 Dyoe 1 Dagy 1 Dye 
D,, ) to model an anisotropic bearing. It has been 


P that rotor elliptical synchronous or 
nonsynchronous response, due to anisotropic bearings, 
is a combination of forward and reverse components, 
which can be measured with two probes. Parameters K', 
D', Ox, Ôp that relate to major/minor stiffness and 
damping axes, cause the reverse component. The 
relationship between the Dynamic Stiffness parameters 
and the bearing coefficients has been established. One 
can easily obtain the bearing coefficients in terms of the 
Dynamic Stiffness parameters, or the Dynamic Stiffness 
parameters in terms of the bearing coefficients. It is 
suggested that the Dynamic Stiffness parameters (which 
give clear insight into instability scenarios), or the 
bearing coefficients as a function of the Dynamic 
Stiffness parameters, be used in rotor dynamic analysis. 
The key parameter, A, which contributes to the 
destabilizing cross stiffness term, A, or Ky, can be 
directly calculated from nondimensional bearing 


coefficients. °R8!T 
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GO ONLINE FOR MORE 
UNDERSTANDING 


The ways in which forward and 
reverse rotating vectors of 
various amplitudes, 
frequencies, and phases 
combine to make circular, 
elliptical, and other orbit 
patterns can be most easily 
understood by using our full 
spectrum simulator, a visual 
tool that allows you to altera 
full spectrum and observe the 
changes in its corresponding 
orbit pattern. 


You'll find the simulator here: 
http://www.bently.com/service! 
training/fss.htm 


learn more online at 





a 


